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ABSTRACT 

The period changes of 86 M5 RR Lyrae stars have been investigated on a one- 
hundred-year time base. The pubhshed observations have been supplemented by 
archival Asiago, Konkoly and Las Campanas photographic observations obtained 
between 1952 and 1993. About two thirds of the O — C diagrams could be fitted 
by a straight line or a parabola. 21 RR Lyrae stars have increasing, 18 decreasing 
and 16 constant period. The mean rates of period change of these variables are: 
(/3) = (P) = -0.006 ± 0.162 dMyr-\ (a) = {p-^P) = -0.021 ± 0.308 Myr'^ Ten 
RR Lyrae stars show fast period decrease with P < — 0.10 dMyr^^. At least some 
of these variables may be in the prc-zero-age horizontal-branch (ZAHB) evolutionary 
stage. The variables on the long-period sequence of the period-amplitude diagram are 
brighter than the other RR Lyrae stars of M5 and are in an advanced evolutionary 
stage moving off from the HB redward. More than one third of the M5 RR Lyrae stars 
investigated have irregular period change. The irregular period behaviour is relatively 
more frequent among the RRc (RRl) stars (50 per cent) than among the RRab (RRO) 
stars (34 per cent). A strict relationship has been found between the irregular period 
change and the Blazhko effect of M5 RRab stars. This fact indicates a common ori- 
gin for these phenomena. It is remarkable that, if the RRab stars showing Blazhko 
effect are omitted from the sample, the mean rates of the period change have small 
positive values ((/3) ^ 0.012 ± 0.147 dMyr^\ (a) = 0.013 ± 0.279 Myr"^), in excellent 
agreement with HB evolutionary model predictions. 

Key words: stars: horizontal branch - stars: oscillations - stars: variables: RR Lyr 
- globular clusters: individual: M5. 



1 INTRODUCTION 

More than ninety years ago lEddingtOT]! (|l918t l had already 
realized that the period changes of a Cepheid-like pulsator 
would give information on the changes of the physics of 
the star's interior during its evolution. The potential of 
measuring the rate of ste llar evolution seemed to be ful- 
filled when iMartinI (| 19381 ) discovered that the periods of 
the RRab stars in omega Centauri were predominantly in- 
creasing. The subsequent observations of RR Lyrae stars in 
other globular clusters, however, failed to demonstrate the 
direct connection between the detected period changes and 
horizontal-branch (HB) evolution. The observed large pe- 
riod decreases, the random and/or abrupt period changes 
could not be reconciled with steady evolutionary effects. It 
has been, however, argued that evolutionary period changes 
must be present in individual stars and must be observable, 
at least on long time-scales. Globular clusters with dozens, 
hundreds of RR Lyrae stars are ideal targets for such stud- 



ies, the mean value of the period-change rates of a large 
sample of RR Lyrae stars of a glob ular cluste r must show 
the trend of their evolution. Indeed. [Led (|l99ll ) proved that 
the mean period-change rates observed in globular clusters 
were consistent with synthetic HB models, if the HB-types 
of the clusters were also considered. 

The generally adopted procedure for investigating the 
period changes of variable stars is the construction and in- 
terpretation of their O — C diagrams. This plots the differ- 
ence between the observed times (O) of a particular phase, 
usually the maximum or the mid-point on the ascending 
branch of the light curve and the predicted time of the same 
phase (C), calculated according to an accepted ephemeris. 
However, using today's computational techniques, the O — C 
data are often calculated as the time shift between the entire 
observed light curve and the normal light curve, rather than 
from one specific phase. This method gives more stable re- 
sults. If the O — C diagram is a linear or quadratic function 
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of time, then the period is constant or continuously changing 
at a constant rate. 

In order to explain the complex nature of some O — C 
diagrams, i.e. the irregular behaviour of the period varia- 
tions that masks the evolutionary changes, several sugges- 
tions have been put forward in the past. The accumulation 
of random variations (some sort of noise) in the period may 
show different drifts in the O — C diagram as a typical real- 
ization of a random walk, even if the period noise has a Gaus- 
sian distribution (see e.g . lSternell934l : iBalazs-Detre fc Petrel 



1 19651 : [Lombard fc Koenj 119931 ) The appearance of this ef- 
fect on the O — C diagram, however, strongly depends on 
the parameters of the probability density function of the 
period noise, and i t does not reflect real period changes. 
ISweigart &: Renzinil (|l979l ) suggested a mechanism for the 
occurrence of irregularities in the pulsation period. Accord- 
ing to their theory, random mixing events through the semi- 
convective zone cause the composition profile of the star's 
core to be altered. This produces changes in its hydrostatic 
structure, resu lting in abrupt a nd continuous period changes 
of both signs. IStotherg (|l980t ) attributed the irregular pe- 
riod changes to hydromagnetic effe cts taking p lace i n the 
atmosphere of the pul sating star. lLaskaride3 (|l974l ') and 
iKoopmann et all (|l994l ) modelled the effect of mass loss on 
the period changes of RR Lyrae stars, while [CouttsI (ligTlal ) 
noticed that binarity might lead to long-term oscillation in 
the star's observed pulsation period. 

Messier 5 (NGC 5904) is one of the globular clusters 
that are richest in RR Lyrae stars. The period changes 
of the cluster's RR Lyrae s tars w e re already investigated 
bv ICoutts fc Sawver Hogp (|l969l). Kukarkin fc Kukarkinal 



l|l97j ). IStorm et all and lReid (19961 ) using baselines 



ranging from seventy to one hundred years. In these previ- 
ous studies, no photometric data were available for nearly 
two decades between 1970 and 1 990, and the most recent 
period-change study (|Reidl Il996l ) concerned only 30 vari- 
ables. Therefore, with the addition of archival data, the re- 
examination of the period changes of a larger sample of vari- 
ables is desirable. 

Now the observational baseline extends over a century, 
and the previously published observations, combined with 
the unpublished ones from the Las Campanas and Konkoly 
Observatories, render almost continuous coverage for the 
last fifty years of the 20th century. This provides the moti- 
vation for the present study of period changes of RR Lyrae 
stars in the globular cluster M5. 



2 DATA AND METHOD OF INVESTIGATION 

2.1 Previous photometric studies 

The first investigation of the variable stars in M5 was carried 
out by iBaUcy (1917|). He made use of 111 plates obtained 
during the years 1895 to 1908 with the 13-inch Boyden, the 
11-inch Draper and the 24-inch Bruce telescopes mostly at 
Arequipa, Peru with exposure times ranging from 20 to 100 
min. This collection was supplemented with 12 photographs 
taken at the Mount Wilson 60-inch refiector with 10-14 min 
exposure times in 1912. Bailey published brightness esti- 
mates for about eighty RR Lyrae stars, and these observa- 
tions provide a firm starting point for the study of period 
changes of the RR Lyrae stars in M5. 



IShaplevI (|l927l ) studied the variables on the basis of 
a collection of 59 plates mostly with double exposures 
(in all 113 exposures) taken with the Mount Wilson 60- 
inch telescope on eight different nights in 1917. Shapley 
did not publish his measurements, and likely they have 
gone astray. Therefore, the plates were measured anew, 
and magnitudes of 62 RR L yrae stars w e re pub lished by 
ICoutts fc Sawver Hoed (|l969l ) and lCouttsI (|l971bh . 

A detailed anal ysis of the variables in M5 was made 
bv lOosterhol (|l94ll ). He investigated the period changes of 
the cluster variables, and noted that the light curve of some 
stars showed RW Dra-type variation (now called Blazhko 
effect). The observations consisted of 74 and 7 plates taken 
with the 60-inch refiector of the Mount Wilson observatory 
in the years 1934 and 1935, respectively. The brightness of 92 
RR Lyrae stars were measured or estimated on this material. 

Comprehensive studies of the period changes 
of RR Lyraes in M5 h ave been carried out by 
ICoutts fc Sawver Hoed (|l969l ) and iKukarkin fc Kukarkinal 
( 197 j ), independentlv of each another. 

Coutts fc Sawyer Hogg published magnitudes of 66 RR 
Lyrae variables of M5 measured on 157 plates taken with the 
74-inch telescope of the David Dunlap Observatory between 
1936 and 1966. 

Kukarkin fc Kukarkina's material (181 plates) was col- 
lected with the 40-cm astrograph at Crimea in the years 
1952-1968. Due to the severe crowding effect near the cen- 
tre of the cluster and the close companions in some cases, 
the brightness of only 51 RR Lyrae variables could be esti- 
mated. 

Since the mid-eighties of the last century, accurate CCD 
observations have been obtained for the RR Lyrae variables 
ofM5^ 

I Cohen fc GordonI (|l987h observed six variables in B 
and / on six nights in 1986. IStorm et al. I ([1991.1 pub- 
lished B and V magnitudes for e leven RR Lyrae star s ob- 
tained during four n i ghts i n 1987. iBrocato et al] (|l996l ) and 
ICohen fc Matthews! ij 19921 ) investigated the M5 variables in 
1989. Cohen fc Matthews made V, I photometry of eight 
stars on six nights, while Brocato et al. gave measurements 
in B, V c o lours for fi fteen variables on four nights. 

iReidI (|l996l ) and iKaluznv et al. I (|2000l ) published com- 
prehensive CCD photometry of the variable stars in M5. 
Reid obtained his V, I observations during ten nights in 1991 
and 1992 and investigated 49 RR Lyrae stars. Kaluzny et al. 
observed 65 RR Lyrae stars of M5 in V colour. Their obser- 
vations were made during 12 nights in 1997. We note here 
that the star identified by Kaluzny as V92 is in fact V17. 

2.2 Observations published in this paper 

In this paper we publish the photographic observations ob- 
tained with the University of Toronto 61-cm telescope at 
the Las Campanas Observatory of the Carnegie Institution 
of Washington between 1972 and 1991, with the 60-cm re- 
flector of the Konkoly Observatory at Budapest during the 
years 1952-1963, with the 1-m RCC telescope of the Konkoly 
Observatory at Piszkesteto Mountain Station from 1976 to 
1993 and with the 1.2-m refiector of the Asiago Observatory 
of University Padova in 1971. 

The Asiago, Konkoly and Las Campanas photographic 
magnitudes of the M5 variables are given in Table [T] avail- 
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able in its full content as Supplementary Material in the 
electronic edition of this article. 

The Las Campanas plates had 103aO emulsion and were 
exposed through a GG 385 filter, so the magnitudes are B 
magnitudes. A total of 508 plates was obtain ed, and they 
were all measured on a Cuffey iris photometer jCuffevlllQSg . 
Il96ll ). The plates from 1972 to 1980 were measured by sum- 
mer students, and the data were used, but not published, 
in the study of the Fourier parameter ip^i bv lClement et all 
1 19921 '). The 1981 to 1991 plates were measured later by C. 
Clement. Since the plates were measured at different times 
by different people, there might be systematic magnitude 
shifts from year to year, particularly for stars that were sub- 
ject to crowding. The 1981-1991 magnitudes should be on 
a consistent system. However, even then, for crowded stars, 
it was not always possible to be consistent because of dif- 
ferences in the seeing conditions. Altogether 73 RR Lyrae 
stars were measured. 

The photographic observations at Budapest with the 
60-cm telescope were exposed on Guilleminot Superfulgur 
or Agfa Astro Spezial plates. Since the sensitivity of the 
emulsions is similar, the measurements or estimates should 
be on a consistent system. On 292 plates the brightness of 39 
RR Lyrae stars was measured. The magnitudes of a further 
35 variables that were close to the crowded central region or 
had close companion(s) were estimated by M. Lovas. 

The observations with the 1-m telescope have been 
made on 103aO plates. In the years 1976, 1988 and 1992-93 
altogether 94 103aO plates were exposed without filters. In 
the other years 104 photographs were obtained on 103aO 
plates through a 2mm GG13 filter. In all, 93 variables were 
measured on these plates. We found no serious systematic 
difference between the filtered and unfiltered blue magni- 
tudes that exceeded the measurement errors. 

The measuring procedure of the photographic plates ob- 
tained with the 60-cm and 1-m telescopes was the follow- 
ing. The plates were digitized on a Umax PowerLook 3000 
flatbed-transparency scanner with 3.13 and 1.29 pixel-per- 
arcsec resolution on the 1-m and 60-cm photographs, re- 
spectively. Digital aperture photometry was applied to the 
images using standard IRAIu packages. In this way, the pho- 
tographic densities of the variables and comparison stars in 
M5 have been determined. 

The brightness of variables on the Asiago plates that 
were obtained in bad seeing conditions have been estimated 
by M. Lovas. These plates had IlaO emulsion and were ex- 
posed without filter. In all 67 variables were estimated on 
35 plates. 

In the course of measuring and estimating the bright- 
ness of the variabl es, the compariso n stars were chosen from 
Arp's photometry (|Arplll955l . ll962l ). His photographic {nipg) 
comparison sequence was applied to the Budapest and Asi- 
ago plates. The Las Campanas and Piszkesteto observations 
are on the B magnitude scale. No systematic difference ex- 
ceeding the 0.02-mag observational error has been found be- 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which arc operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Table 1. Photographic magnitudes of M5 variables. The com- 
plete table is given in the electronic version of the article as sup- 
porting information (Table SI). 

star HJD-2400000 Bpg mag source" 



VI 34122.460 

VI 40985.589 

VI 41447.629 

VI 42598.453 



15.41 
15.80 
15.26 
15.65 



B24 
A48 
L24 
P40 



B24: Budapest 60cm telescope, 
A48: Asiago 1.2m telescope, 
L24: Las Campanas 61cm telescope, 
P40: Piszkesteto Im telescope. 



tween the filtered and unfiltered Piszkesteto blue observa- 
tions. 

These observations supplemented with the published 
ones provide an almost continuous coverage for the second 
half of the last century. 

2.3 The method 

We consider only those RR Lyrae stars of M5 that have a 
long enough historical record to follow their period changes 
i.e. variables with numbering below 100. V53, V86, V93 and 
V94 are omitted because their photometry is seriously af- 
fected by the contamination of close companions. In all, we 
investigate the period changes of 65 RRab and 21 RRc stars 
of M5. 

We used the original photometric data (reviewed in Sec- 
tions 2.1 and 2.2) to derive the O — Cs. Table [5] summarizes 
the available photometry for each RR Lyrae star of M5, the 
sources are denoted by '-I-'. 

The procedure applied to the determination of the O — 
C diagrams and to explore the period c hanges is basically 
the same as used bv lJurcsik et al.l (|200ll ) in the study of cu 
Centauri's variables. 

As the first step, the zero-points of the magnitudes of 
the different observations of each RR Lyrae star were ho- 
mogenized. The mean light curves used for the magnitude 
transformations were the Las Campanas data from 1981- 
1991 and for the 14 stars (V5, V36, V54, V65, V67, V85, 
V87, V90, V91, V95, V96, V97, V98, V99) that were not or 
only partly measured on the Las Campanas plates, they were 
constructed from the 1976-1993 Piszkesteto data. The light 
curves of the RRab and RRc stars have been decomposed 
by fifth and third order Fourier series, respectively, to fix the 
mean light curves. For the sake of uniformity, we used these 
relatively low-order fits for the RRab stars, because higher- 
order solutions gave wavy curves for some of the stars. We 
checked, however, to ensure that the zero-points and phase 
shifts would remain unchanged if higher order fits had been 
used for the variables with good-quality light curves. 

Both vertical and horizontal shifts were determined for 
all the different photometric data sets of the variables in or- 
der to match their normal curves best. The first part of the 
Las Campanas observations, from 1972 to 1980 were treated 
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Table 2. Data. 



Star 


B17 


S27 


041 


Cd9 B24 K71 


A48 L24 P4U C87 d91 


C92 H9d KUU 


VI 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V2 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V3 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V4 


+ 




+ 


+ 


+ + + - - 


+ + 


V5 


+ 




+ 


+ 


+ - + - - 


+ + 


V6 


+ 


+ 


+ 


+ + 


+ + + - - 


+ + 


V7 


+ 


+ 


+ 


+ + + 


+ + + - + 




V8 


+ 


+ 


+ 


+ + + 


+ + + - + 


+ + + 


V9 


+ 


+ 


+ 


+ + + 


+ + + + - 


+ 


VIO 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


Vll 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V12 


+ 


+ 


+ 


+ + + 


+ + + + - 


+ + + 


V13 


+ 




+ 


+ + 


+ + + - - 


+ + 


V14 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V15 


+ 


+ 


+ 


+ + + 


+ + + - + 


+ 


V16 


+ 


+ 


+ 


+ + - 


+ + + - - 


+ + 


V17 


+ 




+ 


+ 


+ + + - - 


+ + 


V18 


+ 


+ 


+ 


+ + + 


+ + + + + 


+ - + 


V19 


+ 


+ 


+ 


+ + + 


+ + + + + 


+ 


V20 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V21 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V24 


+ 




+ 


+ 


+ + + - - 


+ + 


V25 


+ 


+ 


+ 


+ + - 


+ + + - - 




V26 


+ 




+ 


+ 


+ + + - - 




V27 


+ 




+ 


+ + 


+ + + - - 


+ + 


V28 


+ 


+ 


+ 


+ + + 


+ + + + + 


+ + + 


V29 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V30 


+ 


+ 


+ 


+ + + 


+ + + - + 


+ 


V31 


+ 


+ 


+ 


+ + + 


+ + + + + 


+ - + 


V32 


+ 


+ 


+ 


+ + + 


+ + + + + 


+ - + 


V33 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V34 


+ 


+ 


+ 


+ + 


+ + + - - 


+ + 


V35 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V36 


+ 




+ 


+ + 


+ - + - - 


+ 


V37 


+ 




+ 


+ 


+ + + - - 




V38 


+ 


+ 


+ 


+ + 


+ + + - - 


+ + 


V39 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V40 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V41 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V43 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V44 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ + 


V45 


+ 


+ 


+ 


+ + 


+ + + - - 


+ + 


V47 


+ 


+ 


+ 


+ + 


+ + + - - 


+ + 


V52 


+ 


+ 


+ 


+ + 


+ + + - - 


+ + 


V54 


+ 




+ 






+ + 


V55 


+ 


+ 


+ 


+ + + 


+ + + - + 


+ 


V56 


+ 




+ 


+ 


+ + + - - 


+ + 


V57 


+ 




+ 


+ 


+ + + - - 


+ + 


V58 


+ 


+ 


+ 


+ + + 


+ + + - - 




V59 


+ 


+ 


+ 


+ + + 


+ + + + - 


+ + + 


V60 


+ 




+ 


+ 


+ + + - - 




V61 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V62 


+ 


+ 


+ 


+ + + 


+ + + - + 


+ 


V63 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V64 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V65 


+ 


+ 


+ 


+ + + 


+ - + - - 


+ + 


V66 


+ 


+ 


+ 


+ + + 


+ + + - - 


+ 


V67 


+ 


+ 


+ 


+ + + 


+ 




V68 


+ 


+ 


+ 


+ + + 


+ + - - 




V69 


+ 


+ 


+ 


+ + + 


+ + - - 
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Table 2 - continued 



Star B17 S27 041 C69 B24 K71 A48 L24 P40 C87 S91 C92 R96 KOO 



V70 


+ 


+ 4 


+ 


+ 


+ 


V71 


+ 


+ 4 


+ 


+ 


+ 


V72 


+ 


4- 4 


+ 


+ 


4- 


V73 


+ 


4- 4 


+ 


+ 


4- 


V74 


+ 


4- 4 


+ 


+ 


4- 


V75 


+ 


4- 4 


+ 


4- 


4- 


V76 


+ 


4- 4 


4- 


4- 


4- 


V77 


+ 


4- 4 


4- 


4- 


4- 


V78 


+ 


4- 4 


+ 


+ 


+ 


V79 


+ 


+ 4 


+ 


+ 


+ 


V80 


+ 


4- 4 


+ 


+ 




V81 


-(- 


4- -] 


4- 


4- 


4- 


V82 


+ 


4 








V83 


+ 


+ 4 


4- 






V85 


+ 


4 








V87 


+ 


+ 4 


+ 


+ 


4- 


V88 


+ 


4 








V89 


+ 


4 








V90 


+ 


4 








V91 


+ 


4 








V92 


+ 


4- 4 


4- 


+ 




V95 




4 








V96 




4 








V97 




4 








V98 




4 


+ 






V99 




4 










+ 



4- + + - - -4-4- 




References: B17: lBailev| digiTl ): S27: IShaplevI l|l927l ). ICoutts 
ICoutts 



, .. , ■ „ sawyer Hop e* ('1969'), 'Coutts ('l971lj); 041: lOosterhofll l|l94ll '): C69: 

Sawyer Hogd l ll969l ): B24: Budapest 60-cm telescope. TablelT] K71: |kukarkin fc Kukarkina (, igTlj) : A48 : Asiago 1.2 -m telescope 



TablJT] L24: Las Campas 61-cm t e lescop e, Tab ljl] P40: Pis zkeste to 1-m telescope Tabl jl) C87 lCohen fc GordonI lll987l '): S91: IStorm et alj 
il99W . C92: ICohen fc Matthewd l ll992r) ; R96: iReidI lll996l) : KOO: iKaluzny et"aj1 bOOOl') 



as a separate data set because of the zero-point inconsisten- 
cies mentioned in Sect. 2.2. Then, the derived vertical (zero- 
point) shifts have been applied to each data set of each star. 
In this way, consistent data sets have been obtained for the 
variables. The most outlying magnitudes were omitted from 
some of the data sets. 

For homogeneity purposes the B band CCD observa- 
tions were used if they ex i sted. The C CD observations of 
ICohen fc Matthews! (|l992l ). iReidI (|l996l ) and iKaluznv et"al] 
were in V (and /) bands. It is well-known that the 
times of maximum syst e matic ally deviate in blue and yellow 
light, but lJurcsik et all (|20Qll ) have shown if the phase shifts 
of the whole light curve are considered (as in our case), no 
significant difference according to the colours occurs. There- 
fore, CCD V data were used and treated in the same manner 
as photographic and CCD B observations. A similar argu- 
ment holds for the inhomogeneous photographic and CCD 
B data. Although the amplitudes of the light curves and the 
phases of the maximum times might be somewhat difi'erent 
in the different observations, this had a negligible effect on 
the derived phase-shift values that were used for the O — C 
analysis. 

As the second step, the O — C values were determined. 
They have been calculated as the horizontal shifts of differ- 
ent parts of the combined, magnitude zero-point homoge- 
nized data sets to the normal curves. In order to eliminate 
smearing of the light curves caused by period changes, the 



O — Cs should be constructed only from 3-5 years of the 
observations. 

The normal curves that define the zero value of the 
O—C data have been constructed using 100-180 data points. 
Each was checked visually for complete phase coverage. The 
1981-1985 Las Campanas observations were used to con- 
struct the normal curve for most of the variables. For stars 
with rapid period variations (V14, V18, V24, V25, V52, V66, 
V70, V81 and V92) an even shorter time interval (1981- 
1983) has been employed, while data from 1982-1983 de- 
fined the normal curve of V76. For the 14 stars not or only 
partly (V87) measured on the Las Campanas plates, the nor- 
mal curves were constructed from the 1976-1980 Piszkesteto 
data. 

O — C values were calculated, as a rule, for ten groups 
(subsets) of the whole, zero-point corrected data sets of each 
RR Lyrae star. In the case of stars for which measurements 
were not available for all the extended data sources a smaller 
number of groups were composed. Also, if the phase shift 
could only be poorly determined, or the scatter of the light 
curve was too large, groups were merged, in order to stabi- 
lize the solution. Conversely, if the period changed too fast 
and/or irregularly, more than ten groups were formed and 
analysed. The O — C values deduced from the observations of 
iBrocato et al.l (|l996l ) outlie significantly in the O — C plots; 
therefore, they were disregarded in the phase diagrams. 

As starting periods, we accepted the values given in the 
literature (in some cases frequency analysis was needed) and 
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Table 3. Phase shifts and instantaneous periods for VI. The complete table, including all the variables, is given in the electronic version 
of the article as Supporting Information (Table S3). 



VI Pa=0.5217868 d LC JD2444750-2446209" rms = 0.136'' 



time interval [JD— 


Z'-iUUUUUJ 


7 n 

■J u 




w — o [aj 


error (Icr) 




11145.8 - 14437.6 




13570.6 


80 


-0.00084 


0.00275 




14795.6 - 21435.7 




20387.3 


136 


-0.00429 


0.00163 




27540.8 - 31259.6 




28533.9 


128 


-0.00134 


0.00290 




31969.7 - 34929.6 




34124.7 


204 


-0.00449 


0.00114 




35251.5 - 37116.6 




36338.0 


141 


-0.00128 


0.00225 




37436.5 - 39979.4 




38900.4 


193 


-0.00186 


0.00282 




40985.6 - 43286.6 




42383.3 


267 


0.00060 


0.00086 




43631.6 - 45472.5 




44583.2 


287 


-0.00004 


0.00071 




45850.6 - 48013.9 




46684.3 


140 


-0.00244 


0.00099 




48392.8 - 50674.5 




49140.6 


778 


-0.00639 


0.00051 




time interval [JD- 


2400000] 


Td 


N 


P [d] 


error (Icr) 


(P - Pa) * lO^d 


11145.8 - 14437.6 




13570.6 


80 


0.5217821 


.0000023 


-0.47 


14795.6 - 21435.7 




20387.3 


136 


0.5217864 


.0000003 


-0.04 


27540.8 - 31259.6 




28533.9 


128 


0.5217873 


.0000006 


0.05 


31969.7 - 34929.6 




34124.7 


204 


0.5217877 


.0000007 


0.09 


35251.5 - 37116.6 




36338.0 


141 


0.5217860 


.0000023 


-0.08 


37436.5 - 39979.4 




38900.4 


193 


0.5217867 


.0000017 


-0.02 


40985.6 - 43286.6 




42383.3 


267 


0.5217873 


.0000006 


0.05 


43631.6 - 45472.5 




44583.2 


287 


0.5217858 


.0000008 


-0.10 


45850.6 - 48013.9 




46684.3 


140 


0.5217863 


.0000012 


-0.05 


48392.8 - 50674.5 




49140.6 


778 


0.5217858 


.0000002 


-0.10 



" data used to define the normal light curve; LC and PI refer to the Las Campanas and Piszkesteto observations, respectively, 
residual scatter of the time-transformed data with the exclusion of CCD V data. 



constructed preliminary O — C diagrams. Then, the average 
period values (Pa) over the century were determined and 
used in order to obtain a symmetrical form of the diagrams. 
If it was necessary, several trials were done to find the best 
period, which gave the most reasonable O — C plot. This was 
especially the case, if significant abrupt and irregular period 
variations occurred. 

If the O — C is fitted by a polynomial 

k 

0-C = ^c,t'-i (1) 

i = l 

then it can be readily seen that the period variation can be 
eliminated, if the times of the observations are transformed 
according to the equation 

k 

i'=t-Y.c^-\ (2) 

i = l 

The time-transformed data can be coherently phased with 
the period. Pa. 

For some of the variables with strong period varia- 
tions, the zero-point corrections of the different data sets 
had slightly incorrect values because of the smearing of 
the folded hght curves. In such cases, we used the time- 
transformed hght curves to derive refined values for the cor- 
rected magnitudes. 

In summary, the following steps were applied: 

• mo, to — !■ mi, to: magnitude shift of the different obser- 
vations, 

• mi , to — > mi , ti : time transformation of the data using 
equation 2, 



Table 4. The collected photometric observations of M5 variables. 
The complete table is given in the electronic version of the article 
as Supporting Information (Table S4). 



HJD-2400000 


mag 


source" 


magf^ 




star 


11145.794 


14.46 


B17 


14.72 


11145.794 


VI 


21338.885 


14.35 


S27 


14.48 


21338.885 


VI 


11178.785 


14.94 


B17 


15.38 


11178.708 


V2 



" the same as in Table [2] 

* zero-point transformed, homogenized magnitudes 
transformed HJD that eliminates phase/period variations 



• mi,ti — >■ 7712, to: refining the magnitude transformation 
on the time-transformed data, 

• m2,to — >■ 7712, t2: time transformation of the refined 
magnitude data sets, 

where mo, to, mi, ti and 7712, t2 refer to the magnitudes 
and times of the original observations; the zero-point cor- 
rected, time-transformed data; and the refined, zero-point 
corrected, time-transformed data. 

All the obervations used in the analysis and 
the magnitude-zero-point and time-transformed data are 
given in electronic form as Supporting Information 
(Table S4) and can be downloaded from the url: 
http:/ /www. konkoly.hu/24/publications/M5. Table [4] gives 
a sample of the data. 
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O-C [d] 



0.1 




0.1 


-0.1 



(P-Po)*10 
Pn 0.2648176 



O-C [d] 




Po 0.2852373 



Pn 0.2907655 





20000 40000 
JD-2400000 



20000 40000 
JD-2400000 





-0.2 
-0.4 
-0.6 
-0.8 




-0.5 -0.2 0.0 0.2 0.5 
phase 




Pn 0.3213360 




Pn 0.3296029 



Pn 0.3331382 



, ••• • •• 9 



Pq 0.3342683 




20000 40000 
JD-2400000 



20000 40000 
JD-2400000 



-0.5 -0.2 0.0 0.2 
phase 



Figure 1. For each RR Lyrae star, in increasing order of the period, three panels describe the behaviour of the period changes. The left- 
hand panels show the O — C diagrams and their polynomial fits. The O — C values are given in days, and are plotted versus Julian Date. 
In the middle panels the departures from the average values of the period \{P — Pa) X 10^ d] arc plotted from direct period determination 
for the data subsets versus Julian Date. 2a formal errors of the O — C and direct period values arc indicated. For comparison, the 
derivatives of the O — C polynomial fits are also drawn in. The right-hand panels show the folded light curves of the time-transformed 
data calculated from the O — C polynomial fits over the hundred years of observations. The folded light curves give an impression about 
the quality of the O — C fits. 
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Pq 0.4703217 
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Figure 1 - continued 
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Figure 1 - 



continued 
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Figure 1 - continued 
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0-C[d] (P-P„)*10' 0-C[d] (P-P„)*10^ 




20000 40000 20000 40000 -0.5 -0.2 0.0 0.2 0.5 



JD-2400000 JD-2400000 phase 

Figure 1 - continued 



3 O-C DIAGRAMS, PERIOD CHANGES AND 
FOLDED LIGHT CURVES 

The results for the 21 RRc and 65 RRab stars are docu- 
mented in Fig. [T] Three panels are shown for each variable. 
The O — C values and their errors, estimated as twice the 
standard deviation of the corresponding part of the light 
curve from the light-curve template used to determine the 
O — C value, are plotted versus time in the left-side pan- 
els. For the sake of uniformity and simplicity, the O — C 
diagrams have been approximated by polynomials of differ- 
ent order. Whenever it was possible, a linear or quadratic 
fit was applied. Occasionaly, when the period changed very 



irregularly, the O — C diagram was composed from two or 
three polynomials fitted to the different parts of the O — C 
diagram. In these cases, the O — C solution might not be 
unambiguous because of cycle-count uncertainties. 

In the middle panels, the deviations from the average 
values of the periods (P — Pa) x 10^ are plotted from direct 
period determination for the data subsets. 2a formal errors 
of the direct period determinations are also shown. If the 
period could be determined only with large uncertainty, data 
subsets were drawn together. 

The O — C can be given as the function of time as 
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0-C = ^ J P(t)dt - 1. (3) 

Consequently, the temporal period, P(t), can be approxi- 
mated using the Ci coefficients of Eq. 1.: 

k 

p{t) = -p-^^^^T^ +P'^ = Y^i^ - ^w'' + (4) 

For comparison, and to confirm the O — C solution, P{t) 
calculated from the derivatives of the polynomial O — C fits 
are also drawn in the middle panels of Fig [1] 

The right-hand panels show the folded light curves of 
the time-transformed data (Eq. 2.) calculated from the O—C 
polynomial fits over the hundred years of observations. 

Two O — C solutions are shown for some of the stars 
(V8, V19, V27, V38, V62, V65, V78 and V88). The first 
fits the O — C data with linear or parabolic approximation, 
while the second shows a higher order or multiple fit. These 
examples illustrate the strengths of the higher order polyno- 
mial or irregular period variations. We also note here that, 
in some cases (Vll, V17, V56 RRab and V40, V80 RRc 
stars), a sine- like approximation would also be appropriate, 
but for the sake of conformity and uniformity, as mentioned, 
a polynomial fit was always applied. 

Although a visual inspection of the folded light curve 
was already convincing of the validity of both the O — C val- 
ues and their polynomial fits in most cases, we also checked 
this quantitatively. The residuals of the Fourier fits of the 
folded light curves (omitting CCD V data) of the time- 
transformed data have been determined and compared to 
the mean rms of the observations. For variables with period 
variation that was not too complex, the rms of the time- 
transformed data were 0.10-0.15 mag for well resolved stars, 
which is in the same range as the typical errors of the pho- 
tographic observations. The rms of the photographic obser- 
vations of these stars ranged from 0.05 to 0.20-mag (usually 
the smallest residuals were obtained from iris photometry) 
with an average value of about 0.10-0.15 mag. 

The numerical results for the O — C and period data of 
the studied 86 variables are available electronically. Table [3] 
shows an example for VI. The first line of the table identifies 
the variable, gives the period that was used to calculate the 
O — C, the data sets used for constructing the normal curve 
and the rms of the time-transformed light curve. In calcu- 
lating this rms, only photographic and B band observations 
were used. The first part of Table [3] lists the O — C data. 
The time intervals, their mean Julian Date, the number of 
data points, the O — C value and its la error are given. The 
second part of Table [3] gives similar data for direct period 
determinations. The last column of this part of the table 
lists the actual period-deviation values, P — Pa- 

We emphasize here that the Pa periods are not the cur- 
rently best periods of the variables; these periods are the 
best mean periods over the last century. The periods de- 
termined for the latest part of the observations match the 
recent data. 

Amplitude and phase modulations of Blazhko stars 
might introduce some uncertainty into the derived zero- 
point shifts and O — C values in the case of variables 
showing strong light-curve mod ulations. How e ver, m ost of 
the Blazhko stars identified in lJurcsik et al] (|2010l ') show 
strongly irregular period variations, which cannot be used to 



derive any period-change rate. Therefore, these uncertainties 
have no effect on any of the conclusions of this paper. The 
period-change rates of the two Blazhko stars with strong 
phase modulations (V58 and V63), that have dominantly 
continuous period changes are well defined as can be seen 
in Fig. [T] We have to admit, however, that the origin of the 
scatter of the O — C and period variations of these stars may 
either originate from the bias of their phase modulation or 
may reflect real small fluctuations in the periods. 



4 RESULTS AND DISCUSSION 
4.1 Period-change results 

In the previous section the O — C diagrams of 21 RRc and 65 
RRab stars of the cluster are presented. The main charac- 
teristic parameters of the light curves and period changes of 
these stars are collected in Table[Sl In the flrst three columns 
the catalogue numbers of the variable in question, its type 
and the period used to construct the O — C diagram are 
given. 

The next four columns contain some data of the light 
curves that might have a connection with the stars' evo- 
lutionary stage and can be compared with the variables' 
period-change rates: Vi and Vm are the intensity- and 
magnitude-averaged mean V magnitudes, respectively, de- 
rived from the combined R96 and KOO CCD V observations 
for most of the stars. The S91 CCD V phot ometry was used 
for the star V7 and the lBrocato et al.l CCD V data for 
the variables V95 and V98. For Blazhko-stars, light curves of 
the best-represented phase of the modulation were used to 
derive the mean V magnitudes. Bm denotes the magnitude- 
averaged mean B magnitude of the variables determined 
from the CCD B data or, if they are not available then Bm 
has been derived from the Las Campanas or the Piszkesteto 
photographic observations. No mean B magnitude is given 
for V6, V13, V17, V25, V26, V37, V54, V74, V85, V90, V91 
and V92. There are no CCD B observations of these stars, 
and their photographic photometry is seriously affected by 
bright close companions and crowding. The mean photo- 
graphic magnitudes of these stars are about 0.5-1.5 mag 
brighter than expected. The mean-magnitude values were 
calculated as Aq, the flrst member of the Fourier-flt of the 
light curve of the variables. 

The mean V magnitudes of the variables differed by 
0.01-0.03 mag in the R96 and KOO data for the variables 
observed in both surveys. It is thus estimated that the 
mean V magnitudes have similar uncertainties. The error 
of the mean B magnitudes cannot be correctly quantified. 
Although we have omitted variables with strong light con- 
tamination due to crawding, the photographic B magnitudes 
of some of the variables might still be somewhat biased by 
their surroundings. The inhomogeneity of the B data (pho- 
tographic and CCD) introduce further errors. Notwithstand- 
ing these problems, the accuracy of the given B magnitudes 
is estimated to be better than 0.05-mag for most of the stars. 

The pulsation amplitude in the V band. Ay, has been 
derived from the CCD V observations that were used to 
calculate the mean V magnitude. If no V observations 
were available, the Ay amplitude has been estimated as 
Av « Ab/I.Z (|Jurcsik et al.ll2005h . where yls is the blue 
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Table 5. Summary of the light-curve parameters and period-change properties of M5 RR Lyrae stars. 



Star 


type 


-Pa 




* m 








Remarks 


/3 = P 


13 = P 


a = P'^P 


Pmax P^in 






d 












d 


e 


/ 


10-10 


dMyr~i 


lO-iOfj-i 


lO^^d 


VI 


RRab 


0.5217868 


15.103 


15.155 


15.443 


1.07- 


-1.15 


U 


U 


Rl 


0.00 


0.00 


0.00 


0.56 


V2 


RRab 


0.5262647 


15.093 


15.148 


15.559 


1.04- 


-1.11 




9 
Z 


Rl 


- 


- 


- 


1.82 


V3 


RRab 


0.6001852 


15.057 


15.078 


15.451 


0.74 




1 


n 
U 




0.80 


0.029 


1.32 


0.45 


V4 


RRab 


0.4496332 


15.047 


15.093 


15.343 


0.95- 


-1.25 




9 
Z 


Rl 

Dl 


- 


- 


- 


8.22 


V5 


RRab 


0.5458749 


15.107 


15.165 


15.390 


1.02- 


-1.18 


i 


n 
U 


Rl ■ 


-13.11 


-0.479 


-24.02 


5.88 


V6 


RRab 


0.5488297 


15.033 


15.076 


- 


1.00 




i 


U 


Q 


-1.68 


-0.061 


-3.06 


1.08 


V7 


RRab 


0.4943970 


15.035: 


15.103: 


15.491: 


1.30 




1 


U 




7.04 


0.257 


14.23 


2.36 


V8 


RRab 


0.5462315 


15.085 


15.126 


15.497 


0.82- 


-1.08 




1 
i 


Rl 


- 


- 


- 


2.55 


V9 


RRab 


0.6988956 


14.903 


14.931 


15.324 


0.80 




i 


u 




-0.30 


-0.011 


-0.43 


0.17 


VIO 


RRab 


0.5306612 


15.085 


15.133 


15.516 


1.10 




1 


u 




-1.81 


-0.066 


-3.41 


0.61 


Vll 


RRab 


0.5958927 


14.933 


14.988 


15.345 


1.13 






1 

_L 




- 


- 


- 


0.38 


V12 


RRab 


0.4677122 


15.098 


15.172 


15.502 


1.28 






n 




-6.72 


-0.245 


-14.37 


2.44 


V13 


RRab 


0.5131240 


14.955 


15.006 


- 


1.15 




n 
U 


1 
i 
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0.00 


0.00 


0.83 


V14 
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15.424 
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-1.35 




2 


Bl 


- 


- 


- 
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V15 
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15.054 


15.064 


15.355 


0.41 




1 
1 


u 




3.74 


0.137 


11.11 


1.25 


V16 


RRab 


0.6476276 


14.809 


14.869 


15.267 


1.19 






n 
u 




2.06 


0.075 


3.19 


1.03 


V17 


RRab 


0.6013950 


14.942 


14.997 


- 


1.12 




9 


1 
i 


Q 


- 


- 


- 


1.19 


V18 


RRab 


0.4640440 


15.075 


15.151 


15.420 


0.80- 


-1.35 




2 


Bl 


- 


- 


- 
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V19 
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0.4699750 


15.088 


15.173 


15.495 
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9 

z 
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- 


- 


3.84 


V20 


RRab 
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u 
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0.5074854 


- 


- 


- 


- 






9 

z 




- 


- 


- 


18.63 


V26 


RRab 
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0.00 
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15.053 
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- 


- 


- 


1.55 
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0.5439370 
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15.123 


15.511 


0.94 




1 


n 
U 




-6.45 


-0.236 


-11.86 


2.57 


V29 


RRab 


0.4514274 


15.127 


15.164 


15.506 
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-0.88 


i 


i 


Rl ■ 
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-8.01 
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U 
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0.71 
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u 
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Q 


0.00 


0.00 


0.00 
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- 


- 


- 
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Q 
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0.76 
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-5.92 


0.24 


V61 


RRab 


0.5686251 


15.062 


15.095 


15.485 


0.91 




1 







6.44 


0.235 


11.32 


2.41 


V62 


RRc 


0.2814137 


15.082 


15.099 


15.315 


0.53 






2 










1.62 


V63 


RRab 


0.4976794 


15.06: 


15.08: 


15.454 


0.90- 


-1.00 


1 


1 


Bl 


2.22 


0.081 


4.45 


1.25 


V64 


RRab 


0.5445020 


15.093 


15.134 


15.563 


1.02 




1 







-3.79 


-0.138 


-6.96 


1.66 


V65 


RRab 


0.4806696 


15.085 


15.133 


15.415 


0.75- 


-1.00 




2 


Bl 








2.84 


V66 


RRc 


0.3506989 


15.019 


15.031 


15.364 


0.44 






2 










3.53 


V67 


RRc 


0.3490947 






15.137 


0.62 












0.00 


0.00 


0.00 


0.77 


V68 


RRc 


0.3342683 






15.369 


0.50 




1 







-3.96 


-0.145 


-11.84 


1.40 


V69 


RRab 


0.4948735 






15.515 


1.23 




1 







-0.82 


-0.030 


-1.65 


0.49 
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Table 5 — continued 



otar 


type 


PS 


1/6 


1/6 
V 

m 


ab 

ID 

m 






Remarks 


P — F 


fj = F 


OL = 








A 






mag 






d 




J 


10-10 




10-iOd-i 


lO^^d 


V70 


RRab 


0.5585464 


- 


- 


15.576 


1.08 




1 







14.65 


0.535 


26.24 


5.61 


V71 


RRab 


0.5024724 


- 


- 


15.565 


1.23 




1 







4.19 


0.153 


8.33 


1.00 


V72 


RRab 


0.5622688 


- 


- 


15.525 


0.70- 


-0.92 




2 


Bl 


- 


- 




15.11 


V73 


RRc 


0.3401137 


14.958 


14.974 


15.240 


0.51 







1 




0.00 


0.00 


0.00 


0.92 


V74 


RRab 


0.4539926 


15.005 


15.083 


- 


1.39 




1 





d 


-3.40 


-0.124 


-7.49 


1.31 


V75 


RRab 


0.6854167 


14.989 


15.002 


15.498 


0.55 




1 







0.60 


0.022 


0.87 


0.37 


V76 


RRc 


0.4324319 


14.819 


14.827 


15.160 


0.40 






2 




- 


- 


- 


3.63 


V77 


RRab 


0.8451232 


14.760 


14.775 


15.222 


0.60 




1 







5.83 


0.213 


6.90 


3.61 


V78 


RRc 


0.2648176 


15.076 


15.086 


15.332 


0.42 







1 




0.00 


0.00 


0.00 


0.26 


V79 


RRc 


0.3331382 


15.010 


15.020 


15.251 


0.38 












0.00 


0.00 


0.00 


0.33 


V80 


RRc 


0.3365419 


15.042 


15.052 


15.344 


0.41 




2 


1 




- 


- 


- 


0.22 


V81 


RRab 


0.5573072 


15.058 


15.094 


15.524 


0.95 




1 







-13.57 


-0.496 


-24.34 


5.80 


V82 


RRab 


0.5584434 


15.046 


15.079 


15.508 


0.92 




X 







-3.00 


-0.109 


-5.36 


1.19 


V83 


RRab 


0.5533077 


15.092 


15.123 


15.404 


0.84 










d 


0.00 


0.00 


0.00 


0.43 


V85 


RRab 


0.5275226 
















d 








2.83 


V87 


RRab 


0.7383982 


14.915 


14.921 


15.378 


0.37 




1 







6.58 


0.240 


8.91 


4.01 


V88 


RRc 


0.3280850 


15.024 


15.034 


15.305 


0.42 






2 










0.95 


V89 


RRab 


0.5584433 


15.085 


15.122 


15.540 


0.96 




1 







1.26 


0.046 


2.25 


0.61 


V90 


RRab 


0.5571570 












1 





d 


3.14 


0.115 


5.64 


0.99 


V91 


RRab 


0.5849427 


















d 


0.00 


0.00 


0.00 


0.06 


V92 


RRab 


0.4633870 














2 


d 








3.52 


V95 


RRc 


0.2907655 


15.060 


15.079 


15.300 


0.50 
















0.21 


V96 


RRab 


0.5122387 


15.091 


15.136 


15.536 


0.94 








d 








1.69 


V97 


RRab 


0.5446239 


15.042 


15.051 


15.486 


0.40- 


-0.75 






Bl 








1.77 


V98 


RRc 


0.3063745 


15.120 


15.134 


15.401 


0.52 






1 










1.44 


V99 


RRc 


0.3213360 


15.037 


15.049 




0.46 
















0.03 



°': Pa, the period used to construct the O — C. 

^'■Vi, Vm and Bm- intensity and magnitude averaged mean magnitudes. 

Ay. the pulsation amplitude in V. 
Remarks: linear (0), quadratic (1), third order (2), higher order or multiple (-) fit to the O — C. Colon denotes stars with scarce O — C 
data. complexity index of the O — C diagram, straight line or parabola (0), low order fit with scatter or third order fit (1), higher 
order or multiple fit of irregular character (2). This column refers only to RRab stars. Blazhko effect (Bl) and deficient photometry 
for the investigation of light-curve variability (d) are marked. 
3: the full range of the observed period variation. 
See more details in the text. 

amplitude of the light curve derived from the CCD B or pho- 
tographic B observations. For Blazhko stars, the smallest 
and la rgest observed amplitudes, taken from ljurcsik et al.l 
l|2010l . Paper II), are given. 

As the inhomogeneity of the listed mean magnitudes 
and amplitudes gives rise to enhanced uncertainties, these 
values have to be taken with caution and are given only for 
guidance. As no accurate multicolour photometry of the full 
sample of M5 variables exists, this combined photometric 
information is, however, also utilized on a statistical ground. 

The 'Remarks' supply information on the complexity 
and quality of the O — C diagram, i.e. on the derived period 
changes and their reliability for each variable, and on the 
presence of light-curve variation for RRab stars. The first 
column of the 'Remarks' gives the order of the polynomial 
fit to the O — C data. The index '0' means straight-line fit, 
'1' means quadratic, while '2' means cubic approximation. In 
this column, a 'hyphen' indicates that an even higher-order 
polynomial is fitted (e.g. in case of sine-like O — C variations) 
or a combined polynomial approximation is applied to sat- 
isfactorily describe irregular period variations with sudden 
changes. In these latter cases, when the O — C cannot be 
fitted with a single polynomial, cycle-count uncertainties of 
the O — C solution cannot be excluded. As these O — Cs 



are always classified as irregular, i.e. no period-change rate 
is derived from these data, these ambiguities have no effect 
on any of our conclusions. 

There are five stars (V85, V95, V96, V97 and V98) de- 
noted by colons in this column that refer to their scarcely 
defined O — C diagrams. These stars are included only for 
completeness, but the O — C solutions shown in Fig. [T] are 
ill-defined, other solutions are similarly possible. These stars 
are disregarded in the following. 

The second column of the 'Remarks' gives information 
about the complexity of the O — C diagram, i.e. about the 
irregularity of the period changes. The irregularity indices '1' 
and '2' refer to moderate and strong irregularity of the O — C 
diagram while '0' means that the variation in the pulsation 
period is smooth and regular, i.e. it is either constant or 
linearly changing. Sometimes the distinction between the 
classes '0' or '1' and '1' or '2' may seem arbitrary. This, 
however, does not affect the subsequent discussion. 

The third column of the 'Remarks' refers to the pres- 
ence of light-curve variation for fundamental-mode RRab 
stars (see details in Paper II). 'Bl' denotes that the star 
definitely exhibits Blazhko effect, while 'Bl:' indicates that 
light-curve modulation is probable. RRab stars whose pho- 
tometry is seriously affected by crowding or close, bright, 
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neighbouring star(s) are labelled by 'd'. The deficiency of 
the photometry of these stars makes any statement on their 
light-curve variation impossible. 

The last four columns of Table [5] give information on 
the observed period changes of the cluster's RR Lyrae stars: 
/3 = AP/At in 10"^" and dMyr"^ and a = P'^AP/At in 
10^^"d^^ units, while Pmax — Pmin is the full range (the 
maximum fluctuation) of the observed period variation. 

We have to admit that in some cases it is difficult to 
decide whether the period is constant with some random 
noise or it has a very slow monotonic change. If the absolute 
value of the coefficient of the quadratic term in the parabolic 
approximation was smaller than its 2a error then the O — 
C diagram was fitted by a straight line, and /3 = was 
accepted. Of course, for these stars the formally calculated 
full range of the period variation may differ slightly from 
zero due to small scatter originating from real period noise 
or observational error. 



4.2 Evolutionary considerations 

The final goal of this study is to determine if the observed pe- 
riod changes of the M5 RR Lyrae stars could be attributed, 
at least in the mean, to stellar evolution. From evolution- 
ary consideration it is well-known that the period increases 
as a star evolves redward, and decreases as it moves blue- 
ward in the HR diagram. Theory predicts evolution in both 
directions and with different period-change rates depend- 
ing on the HB-evolutionary stage of the RR Lyrae stars. 
It has also been known for some time that RR Lyrae pe- 
riod changes may be characterized by strong irregularities. 
As Table [5] shows, more than one third of the O — C dia- 
grams of M5 RR Lyrae stars cannot be interpreted solely 
by constant or linearly changing periods. In spite of the fact 
that evolutionary changes are masked by irregularities in a 
number of cases, it is hoped that the period-change rates 
of RR Lyrae stars with steady period variation give some 
information on their evolution. 

The mean and the most-probable values of the period- 
change rates for the RRc, RRab and the total sample of RR 
Lyrae stars that have straight-line or parabolic O — C dia- 
grams are summarized in the first part of Table [6l Although 
the averages of the a and /3 values suggest a slight, overall 
period decrease, it is not the least significant. The most- 
probable values of the period changes of the statistically 
larger samples (RRab and all the variables) are practically 
zero. 

The generally accepted HB-evolutionary models do not 
indicate that high, negative period-change rates are likely, 
and large, positive values are supposed to occur only in the 
very late stages of HE evolution at luminosities brighter 
than on the zero-age horizontal branch (ZAHB). The period- 
change rate oi versus period for the M5 variables is shown 
in Fig. [21 and the largest value of the normalised period 
variation versus period is plotted in Fig. [S] From these 
figures it can be immediately seen that in M5 (like in 
other globular clusters) there are variables with extreme a 
values. These extreme period changes occur among RRab 
rather than among RR c stars according to the study of 
iRathbun fc SmithI |l997l ). 

In our sample, there are three RRab stars with extreme 
period-change rates, two stars (V5 and V81) with large neg- 



Table 6. Statistics of the pulsation-period-change rates of RR 
Lyrae stars in M5 





/3[dMyr-l] 


a [Myr^i] 


No. 


total sample 


RRc mean 


-0.003 ± 0.073 


-0.014 + 0.221 


11 


m.p.v." 


-0.006 


-0.021 




RRab mean 


-0.006 ±0.178 


-0.023 + 0.328 


44 


m.p.v. 


4-0.003 


+0.001 




all mean 


-0.006 ±0.162 


-0.021 + 0.308 


55 


m.p.v 


+0.002 


-0.002 




without 3 outliers 


RRab mean 


+0.004 ±0.123 


-0.005 + 0.234 


41 


m.p.v. 


+0.008 


+0.012 




all mean 


+0.003 + 0.114 


-0.007 + 0.229 


52 


m.p.v 


+0.005 


+0.005 




without Blazhko RRab stars 


RRab mean 


+0.017 + 0.163 


+0.021 + 0.296 


38 


m.p.v. 


+0.015 


+0.024 




all mean 


+0.012 + 0.147 


+0.013 + 0.279 


49 


m.p.v 


+0.009 


+0.014 





" the most-probable value, i.e. the mode 
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Figure 2. Period-change rates ( a = P/P[10~^''d-i] ) of the RR 
Lyrae stars in M5. Circles and asterisks denote fundamental-mode 
and first-overtone variables, respectively. RRab stars showing the 
Blazhko effect are shown by filled symbols. Variables with irreg- 
ular period-change behaviour are set at o = —30 to show their 
period distribution. Positive and negative period-change rates are 
equally frequent among M5 variables, with extreme (|q!| > 20) 
values for throe RRab stars. Period decreases, irregular period 
change and Blazhko modulation occur in RRab stars only with 
period shorter than 0.6 d. 
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Figure 3. Maximum values of the normalised period variation 
of the RR Lyrae stars in M5. Filled circles and diamonds denote 
variables with strong and moderate irregularity of their period 
changes, stars with regular period variations (constant or lin- 
ear) are shown by open circles, asterisks are for variables with 
not enough information on the period-change properties. Period 
changes of the order of 10~* d are observed in both RRab and 
RRc stars. These very large values of period changes are detected 
only in variables with irregular period changes and in RRab star 
only at shorter periods. 

ative and one (V70) with large positive values. V70 lies far 
from the cluster centre, so it was not included in any of 
the CCD studies. Therefore, only its photographic B mag- 
nitude is available. The mean Bpg brightness of V70 is one 
of the faintest among the whole sample, indicating that V70 
cannot be in a highly evolved evolutionary phase. We thus 
suppose that the period changes of these three stars are not 
induced by evolutionary effects. Since their extreme values 
strongly influence the mean period-change rates, the aver- 
ages and the most-probable values of the period changes are 
also derived excluding them (see the second part of Table [6]). 
The most-probable a and /3 values for RRab and for all the 
52 stars correspond to small positive rates, while the mean 
values have very small positive and negative /3 and a values, 
respectively. 

We thus conclude that, on the average, the period- 
change rate of RR Lyrae stars in M5 are very close to zero. 
We also note that the median rates of the period changes for 
all the different samples is zero . This conclusion confirms the 
result of th e previous studies jCoutts fc Sawyer Hogglll969l : 
lReidlll996l ) on an extended baseline and larger sample, and 
is in good agreement with the prediction of Lee's (1991) cal- 
culation on the stellar evolution - period change connection: 
the mean (median) rate of period change should be close to 
zero (0.005-0.010 dMyr~^) for clusters with HB-type like 
M5. 

An other aspect of the problem of period changes is the 
distribution of the direction of the changes, whether the pro- 
portion of the decreasing and increasing periods reflects the 
evolutionary state of the RR Lyrae stars of a globular clus- 
ter. The mod els of post-ZAHB evolution (|Lee et all Il990l : 
iDormarJl 19921 ) predict that after an initial 'hook' to the red 
caused by shell readjustment, RR Lyrae stars evolve blue- 
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beta[10"'°l alpha[10"'°d"'] 

Figure 4. Histograms of the distributions of the period-change 
rates /3 [lO-^"] and a [lO'^^A-'^]. RRab and RRc stars are shown 
by gray and black areas, respectively. Although there is some hint 
that small positive period values are in excess to negative ones 
in RRab stars, this effect statistically is only weakly significant. 
On the average, the period-change rates of M5 variables are very 
close to zero. 



ward (with decreasing periods) across the instability strip, 
and that when they have depleted the helium in their cores, 
they increase in luminosity and evolve redward (with in- 
creasing periods). The histograms of the distribution of the 
period-change rates jS and a are shown in Fig. |4] This fig- 
ure suggests a slight excess of RRab stars with small in- 
creasing period. Our sample has 19 ab stars with increasing 
period, contrasting with 14 ab variables with decreasing pe- 
riod, while 11 ab stars have constant period. Among the 
RRc stars 2 have increasing, 4 decreasing and 5 constant 
period. Hence, the full sample contains 21 variables with in- 
creasing, 18 with decreasing and 16 with constant period. 
The small number of the c-type stars obviously precludes 
the possibility of their distinct statistical treatment. 

A simple test can prove that the observed small excess of 
RRab stars with increasing period is only weakly significant. 
If we distribute equally the stars with /3 = 0, we end up 
with 24 or 25 positive and 20 or 19 negative. Supposing that 
the likelihood of the positive and negative period change is 
equal, the expected value of the binomial distribution is 22 
with a variance of 3.32. For the whole sample, the expected 
value is 27.5 with a variance of 3.71, so the distribution of 
the P value signs does not provide clear statistical evidence 
for an excess of period increases among all the M5 RR Lyrae 
stars. This might be a consequence of the fact that the mean 
period-change rate of the variables is in agreement with the 

~zero mean period-change rates of the v ariables. 

In M3, another Ool-type cluster, ICorwin fc CarnevI 

(|200lf ) found 21 RRab stars with increasing and 14 with 
decreasing period, i.e. a small excess of increasing periods. 
This might be explained by that a larger fr action of the 
RRab stars of M3 shows OoII characteristics (|jurcsik et al.l 
l2003l ) than that of the M5 RRab variables. At the same 
time, the RR Lyr ae stars of the well-observ ed OoII-type 
clusters like M15 (ISilbermann fc SmithI [1995I ) and uj Cen 
(|jurcsik et al.l [200J) show a definite surplus of increas- 
ing periods. This implies that the majority of the vari- 
ables in these clusters are evolving across the instability 
strip from blue to red. All these results strengthen that 
the Oosterhoff properties of globular-cluster variables are 
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Figure 5. Intensity-averaged mean V briglitness versus period 
diagram of tlie RR Lyrae stars in M5. Symbols are the same as 
in Fig. |3] The period-change rate and its direction are indicated 
by horizontal lines. RRab stars with period decrease are char- 
acterized by short period and fainter luminosities, large period 
increase occurs at longer period, more luminous variables. 



stronR ly connected to their hor izontal-branch evolutionary 
state JClement fc Sheltonlll999|V_ Nevertheless, it should be 
noted that ISmith fc Sandagel (|l98lh called attention to the 



fact that no cluster exhibits a large plurality of decreasing 
period. 

On the whole, the HB evolution is a plausible expla- 
nation for increasing periods. Stars within the instability 
strip evolve rapidly from blue to red, toward the end of core 
helium burning, producing a large positive rate of period 
change. In this respect, it is worth mentioning the difference 
between the period-change behaviour of the RRab stars of 
the clusters M3 and M5. Although both clusters belong to 
Oosterhoff-type I and resemble each other in many aspects, 
M5 has fewer strongly increasing periods than M3. In M5, 
only five out of the 44 RRab (11 per cent) have very large 
period-increase rate (/3 > 0.20 dMyr"^), while in M3, 29 
per ce nt of the RRab stars have . In oth er Ool-type clusters , 
M14 dWehlau fc FroelichI Il994l). M 28 (|Wehlau et al.1 Il986l ) 
and NGC 7006 (|Wehlau et al.lll999l ) the frequency of RRab 
stars with strongly increasing periods are 16, 8 and 10 per 
cent, respectively. This simple statistic shows that the M5 
RRab stars behave Uke other Ool-type clusters, and only 
MS behaves somewhat differently. 

One of the most intriguing issues connected with the 
period changes of RR Lyrae stars is the large negative val- 
ues of /3, observed in different clusters. Among the 44 RRab 
stars of M5 that have parabolic or straight-line O — C dia- 
grams, nine (20 per cent) show decreasing periods with as 
high rate as /3 < —0.1 dMyr~^, that certainly cannot be 
explained by canonical HB evolution. It is worth mention- 
ing that the mean period of the sample of 44 RRab stars is 
0.5622 d, whereas the mean period of the nine ab stars with 
fast decreasing period is significantly less, 0.5127 d. 

A comparison with other well-studied globular clus- 
ters proves that the occurrence of RRab stars with very 
fast period decreases is a general feature, characteristic 



of both Oosterhoff-types. For ex ample, 11 out of the 35 
RRab stars (31 per cent) in M3 jCorwin fc Carnevl [200lh 
and 9 out of the 42 ab va riables (21 per cent) in NGC 
6934 jStagg fc Wehlaulll980l l exhibit large period decreases 
{P < —0.10 dMyr"^). The OoII-type clusters also have (per- 
haps relatively less in numbers) RRab stars with strong pe- 
riod decreases. In uj Cen, 29 RRab stars of the chemically 
homogeneous group have linear period changes, and among 
the m, four (14 per cent) have strong period decreases (table 
7 m ljurcsik et alll200ll ). In M15, there are two out of 13 
RR ab stars (15 per cent) th at show strong period decreases 
( Sil bermann fc SmithI [l995l ) . The trend that the mean pe- 
riod of the ab stars with fast period decreases is shorter 
than the mean period of the whole sample also holds for the 
clusters mentioned. 

To resolve the contradition between canonical HB evo- 
lutionary model predictions and the observed frequency of 
strong period decreases, it is generally supposed that most 
of the period changes are d ue to some kind of 'noise' r ather 
than to evolutionary effects. [Sweigart fc Renzinil (| 19791 ) have 
found that the mixing events at the convective core edge, 
the transfer of helium into the convective core and the sub- 
sequent chemical readjustment of the semi-con vective zone 
around the core, is an intrinsically noisy process that can 
lead to large period-decrease rates . 

Recently, ISilva Aguirre et al.l (|2008l ) have given an al- 
ternative explanation for the frequent occurrence of strong 
period decreases among the cluster variables. They con- 
structed pre-ZAHB evolutionary tracks for a chemical com- 
position appropriate to the globular cluster M3, and inves- 
tigated the period-change behaviour of variables in the final 
approach to the ZAHB location. They have found that, be- 
fore settling on the ZAHB, the variables are subject to a 
strong period decrease with the most likely P values around 
—0.3 dMyr^^, but more extreme values (—0.8 dMyr~^) may 
also take place. The model simulations have also shown that 
some percent of the RR Lyrae population are in the pre- 
ZAHB evolutionary state, and the pre-ZAHB pulsators are 
expected to have longer periods than the bona fide HB pul- 
sators. As M5 resembles M3 in different asp e cts, t herefore 
the model calculations of ISilva Aguirre et al.l (|2008h can be 
compared with our results. The high percentage of variables 
with decreasing periods and the relatively shorter mean pe- 
riods of these stars in M5 show that, likely, the pre-ZAHB 
evolution cannot fully explain the observations. Neverthe- 
less, further model calculations of the pre-ZAHB evolution- 
ary phase may lead to more satisfactory results. 

It is perplexing, however, why the observed period- 
increase rates are in very good agreement with evolutionary 
predictions, if the period-decrease rates are not. If this is in- 
deed the case, some significant difference between RR Lyrae 
stars with increasing and decreasing period rates should be 
found, but observationally, there is no evidence of any dif- 
ference between these stars. 

It is an interesting question whether the evolutionary 
effects are apparent in the different diagrams of cluster 
variables, e.g. period-amplitude, period-brightness, etc. di- 
agrams. Fig. [5] and Fig. |6] show the relationships between 
the intensity-averaged mean V brightness {V) and period, 
as well as between the V amplitude {Ay) and the period. 
The period-change rate and its direction are indicated by 
horizontal lines in both figures. The (V) versus period plot 
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Figure 6. V amplitude versus period plots of the RR Lyrae stars in M5. In the left-hand panel symbols indicate the regular and 
irregular period-change behaviour of the stars as in Fig. [5] In the right-hand panel the symbols indicate the intensity-averaged mean V 
brightness of the variables, open circle: (V) < 15™0; filled circle: 15™0 < {V) < 15™ 1; filled diamond: {V) > 15™ 1. The period-change 
rate and its direction is shown by horizontal lines in both panels. In the left-hand panel we can also recognize that the period changes of 
short-period, large-amplitude RRab stars are often irregular. Large negative period changes occur in short-period RRab stars at fainter 
mean magnitudes while large period increases arc observed at the whole period range of RRab stars, but the amplitudes of these variables 
tend to be larger than that of the majority of the stars, indicating their evolved HB stage. 



displays that the variables with longer periods have higher 
luminosity and constant or increasing period (the only ex- 
ception is V9 with a very slow = —0.011 dMyr"^ de- 
crease). The long-period sequence in the Ay versus period 
diagram shows the same trend. In accordance with evolu- 
tionary theories these stars are at an advanced e volutionary 
stage moving off from the HB (|Clement fc Shelton 1999 ). 

One of the most important issues that should be ad- 
dressed, is the irregular period changes of RR Lyrae stars 
in globular clusters. Scrutinizing the O — C diagrams of 81 
RR Lyrae stars in M5, it turned out that 31 out of them (38 
per cent) cannot be satisfactorily fitted by straight line or 
parabola. According to the number of inflection points and 
the order of the polynomial fit, these diagrams were classi- 
fied as moderately or strongly irregular (see index '1' and 
'2' in the middle column of the 'Remarks' in Table [Sjl. If 
we separately investigate the irregular period behaviour of 
the fundamental-mode and first-overtone pulsators, we find 
that more than the half of the c-type stars (10/18) have 
complex O — C diagrams, while one third of the RRab stars 
(21/63) have irregular period changes. The fact that irregu- 
lar period behaviour is more frequent a mong RRc sta r s than 
RRab stars has already been noted bv lJurcsik et al.l (|200ll ) 
in a similar study of the period-change behaviour of the uj 
Cen variables. A comparison with other clusters' RRab and 
RRc stars prove that thi s period behaviour tends to be a 
general feat ure. In M14 (IWehlau &l Froelichl [19941 ) and in 
NGC 7006 (|Wehlau et al.l ll999l ) the frequency of irregular 
period changes of ab-stars are 20 and 14 percentages, while 
those of the c-type stars are 80 and 43 percentages, respec- 
tively. 

The most generally accepted explanation for the irreg- 



ular p eriod changes was put forward bv lSweigart fc Renzinil 
The mixing process within the semi-convective zone 
and the overshooting at the convective-core edge were found 
to be responsible for the non-evolutionary changes. These 
physical processes, however, cannot account for the different 
period-change behaviour of the RRab and RRc stars. An- 
other problem is that the observed irregular period changes 
are as large as 2-5 -10"* APP^^ in some cases (e.g. in V4, 
V14, V18, V52 and V72) occurring on time-scales of some 
tens of years only. Although many parameters influence the 
time-scales and the resultant period-change values caused 
by mixing events, as an average value, 370 yr was derived 
for the time interval required to produce a period change o f 
3 • 10"^ (APP~^) (see flg. 9 in Isweigart fc Renzinil [l979l ) . 
The two orders of magnitude discrepancy between observa- 
tions and model predictions may give rise to some doubt 
that the rapid, irregular period variations of RR Lyrae stars 
can indeed be generated by the proposed mixing events in 
the interior. To settle this issue, more detailed numerical 
modelling of the effect of the mixing events on the pulsation 
period is highly needed. 

Here we outline another approach to the interpretation 
of irregular period changes of RR Lyrae stars. Fig. [7] shows 
that there is a striking connection between the irregular pe- 
riod changes and the Blazhko effect of M5 RR Lyrae stars. 
We omit the RRc stars from the following discussion, the 
detection of Blazhko effect in these stars is difficult as the 
modulation amplitude is usually in the range of the pho- 
tometric accuracy. The figure demonstrates that all RRab 
stars exhibiting very complex O — C diagram, i.e. irregular 
period changes, have been found to show the Blazhko effect, 
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Figure 7. Connection of irregular period changes with Blazhko 
effect of RRab stars. Stars with strong and moderate irregular pe- 
riod changes are set at 2 and 1 values. Stars with constant period 
and with linear period changes are plotted at 0. Filled and open 
circles denote Blazhko and mono-periodic RRab stars. Variables, 
which have scarce and defective photometry arc denoted by as- 
terisks. Light-curve modulation of these stars cannot be detected 
because of the inaccuracy of their photometry. There is a striking 
coincidence of variables with complex period-change behaviour 
and with the occurrence of light-curve modulation. 



if the sample of variables with good enough photometry to 
detect light-curve variability is considered. 

It seems that the inverse statement is also true with only 
a few exceptions. Two Blazhko stars, VI and V30 have con- 
stant period, and two, V5 and V58 have decreasing period 
with very high rates, -0.479 dMyr"^ and -0.219 dMyr"\ 
respectively. If, however, the O — C diagrams of field Blazhko 
RR Lyrae stars are considered, it can be immediately real- 
ized that the irregular (from time to time abrupt) period 
changes are followed by a long interval of linearly changing 
period, a period readjustment (see e. g. the O — C diagrams o f 
XZ Cyg, RR Gem, SZ Hya in fig. 4 of lLe Borgne et al.ll200'if ). 
On the basis of all this, the period changes of Blazhko RR 
Lyrae stars cannot be identified with evolutionary effects, 
and should be left out of evolutionary considerations. 

If variables showing Blazhko effect are omitted then the 
period-change rates of RRab stars with constant or linearly 
changing period (38 stars) have small positive mean and 
most-probable values. If RRc stars are added to this sam- 
ple (49 stars) the mean and the most-probable values of 
the period-change rates hardly change. The last section of 
Table |B] summarizes the mean period-change rates and the 
most-probable values if the Blazhko RRab stars are left out 
of consideration. Curiously enough, these value s are in ex - 
cellent agreement with HB-model predictions of lLed (|l99lf ). 

The distribution of the period-change rates of the RRab 
stars also changes when the Blazhko stars are omitted. Now, 
among the 38 RRab stars 11, 9 and 18 have decreasing, 
constant and increasing period-change rates, respectively. 
In this case, a binomial distribution with equal likelihood of 
both period-change directions has a variance of 3.1, while 
the difference between period changes is 7, more than twice 



as large. So, in this sample, the period of RRab stars is 
dominantly increasing. 



5 SUMMARY 

The period changes of eighty-six M5 RR Lyrae stars have 
been investigated. The published data have been supple- 
mented by archival, previously unpublished observations 
covering the second half of the last century. The observa- 
tions of M5 variables now provide an almost continuous time 
coverage for about a hundred-year time base, making the ac- 
curate study of the period changes possible. 

Reliable O — C diagrams (phase diagrams) were con- 
structed for 81 RR Lyrae, 62 RRab and 19 RRc stars. The 
O — C plots of 44 RRab and 11 RRc stars could be well 
approximated by a straight line or parabola, indicating con- 
stancy or linear change of the pulsation period of these stars 
on the hundred-year time base of the observations. 21 RR 
Lyrae have increasing, 18 decreasing and 16 constant pe- 
riod. The excess of increasing period among RRab stars is 
slightly more pronounced, however, statistically that is still 
weakly significant. Omitting Blazhko variables the situation 
changes, there are 18 stars with increasing, 9 with constant 
and only 11 with decreasing period in this sample. 

The average, median and most-probable values of the 
period-change rates of M5 variables equals zero, or have very 
small, either positive or negative values for different sub- 
samples. Omitting the Blazhko stars, however, the period- 
change rates indicate a small period increase on the average. 
Its most-probable value for 49 stars is 0.009 dMyr^^, in very 
good accordance with predictions of synthetic HB-model re- 
sults (,LcQil991l) . 

Large period-change rates (both negative and positive) 
characterize some of the shorter-period RRab stars, while 
variables with period longer than 0.6 d always have constant 
or increasing period. The observed rates of period increase 
are, in general, in good agreement with standard HB evolu- 
tionary model predictions, but most of the period-decrease 
rates are significantly larger than tha t canonical models 
would allow. The recent supposition of ISilva Aguirre et al.l 
(|2008l ) that connects these changes to pre-ZAHB evolution 
is a promising attempt, but fails to explain the statistics and 
the period distribution of these stars. 

More than one third of the O — C diagrams of M5 
RR Lyrae stars cannot be fitted solely by straight line or 
parabola, i.e. the period changes of these stars are moder- 
ately or strongly irregular. The pulsation periods of some 
of these stars vary by some 10"'* d within relatively short 
time intervals (10-30 yr). The irregular behaviour of the pe- 
riods is more frequent among RRc stars, and it is rather 
characteristic of the shorter-period RRab stars. An impor- 
tant relation between the irregular period change and the 
Blazhko effect has been revealed: in M5, stars with strongly 
irregular period changes always have Blazhko modulation. 



ACKNOWLEDGMENTS 

The constructive, helpful comments of the referee, Katrien 
Kolenberg are much appreciated. We thank Dr. Katalin 



20 Szeidl et al. 



Olah for placing her photographic observations at our dis- 
posal. The financial support of OTKA grant K-068626 is 
acknowledged. C. Clement thanks the Natural Sciences and 
Engineering Research Council of Canada for financial sup- 
port. Zs. H. thanks the 'Lendiilet' program of the Hungarian 
Academy of Sciences for supporting his work. 



Storm, J., Carney, B.W., Beck, J.A., 1991, PASP, 103, 1264 
Stothers, R., 1980, PASP, 92, 475 
Sweigart, A.V., Renzini, A., 1979, A&A, 71, 66 
Wehlau, A. et al., 1986, AJ, 91, 1340 
Wehlau, A., Froelich, N., 1994, AJ, 108, 134 
Wehlau, A., Slawson, R.W., Nemec, J. M., 1999, AJ, 108, 
134 



REFERENCES 

Arp, H., 1955, AJ, 60, 317 

Arp, H., 1962, ApJ, 135, 311 

Bailey, S.I., 1917, Harvard Ann., 78, 103 

Balazs-Detre, J., Detre, L., 1965, in Strohmeier, W. ed. 

The Position of Variable Stars in the Hertzsprung-Russell 

diagram, Proc. lAU Coll. No. 3, Veroff. der Remeis- 

Sternwarte Bamberg IV. (Nr. 40), 184 
Brocato, E., Castellani, V., Ripepi, V., 1996, AJ, 111, 809 
Clement, C, Jankulak, M., Simon, N. R., 1992, ApJ, 395, 

192 

Clement, C, Shelton, I., 1999, ApJ, 515, L85 
Cohen, J.G., Gordon, G.A., 1987, ApJ, 318, 215 
Cohen, J.G., Matthews, K., 1992, PASP, 104, 1205 
Corwin, T. M., Carney B. W., 2001, AJ, 122, 3138 
Coutts, Ch. M., 1971a, in Strohmeier, W. ed. New Direc- 
tions and New Frontiers in Variable Star Research, Proc. 
lAU Coll. No. 15, Veroff. der Remeis-Sternwarte Bamberg 
IX.(Nr.lOO.), 238 
Coutts, Ch. M., 1971b, Publ. David Dunlap Observatory 
3(3), 79 

Coutts, Ch. M., Sawyer Hogg, H., 1969, Publ. David Dun- 
lap Observatory 3(1), 1 
Cuffey, J., 1956, Sky and Telescope, 15, 258 
Cuffey, J., 1961, AJ, 66, 71 
Dorman, B., 1992, ApJS, 81, 221 
Eddington, A. S., 1918, MNRAS, 79, 2 
Jurcsik, J., Sodor, A., Varadi, M., 2005, IBVS, No.5666 
Jurcsik, J. et al., 2001, AJ, 121, 951 
Jurcsik, J. et al., 2003, ApJ, L597 

Jurcsik, J. et al., 2010, MNRAS, submitted (Paper II) 

Kaluzny, J. et al., 2000, A&AS, 143, 215 

Koopmann, R., Lee, Y.-W., Demarque, P., Howard, J. M., 

1994, ApJ, 423, 380 
Kukarkin, B.V., Kukarkina, N.P., 1971, Perem. Zvezdy 

Pril., l(No.l), 1 
Laskarides, P. G., 1974, Ap&SS, 27, 485 
Le Borgne, J.-F. L. et al., 2007, A&A, 476, 307 
Lee, Y.-W., 1991, ApJ, 367, 524 
Lee, Y.-W., Demarque, P., 1990, ApJS, 73, 709 
Lee, Y.-W., Demarque, P., Zinn, R., 1990, ApJ, 350, 155 
Lombard, F., Koen, Ch., 1993, MNRAS, 263, 309 
Martin, W. Chr., 1938, Leiden Ann., 17, No.l 
Oosterhoff, P. Th., 1941, Leiden Ann., 17, No.4 
Rathbun, P. G., Smith, H. A., 1997, PASP, 109, 1128 
Reid, I.N., 1996, MNRAS, 278, 367, Erratum: 1996 MN- 
RAS, 282, 304 
Shapley, H., 1927, Harvard Bull., 851, 15 
Silbermann, N. A., Smith, H. A., 1995, AJ, 109, 1119 
Silva Aguirre, V. et al., 2008, A&A, 489, 1201 
Smith, H. A., Sandage, A., 1981, AJ, 86, 1870 
Stagg, Ch., Wehlau, A., 1980, AJ, 85, 1182 
Sterne, Th.E., 1934, Harvard Obs. Circ, No.386-387 



